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ABSTRACT 

We solve for the local vertical structure of a thin accretion disk threaded 
by a poloidal magnetic field. The angular velocity deviates from the Keplerian 
value as a result of the radial Lorentz force, but is constant on magnetic 
surfaces. Angular momentum transport and energy dissipation in the disk 
are parametrized by an a-prescription, and a Kramers opacity law is assumed 
to hold. We also determine the stabihty of the equilibria with respect to 
the magnetorotational (or Balbus-Hawley) instability. If the magnetic field 
is sufficiently strong, stable equilibria can be found in which the angle of 
inclination, i, of the magnetic field to the vertical at the surface of the disk has 
any value in the range < i < 90°. By analyzing the dynamics of a transonic 
outflow in the corona of the disk, we show that a certain potential difference 
must be overcome even when i > 30°. We determine this potential difference as 
a function of i for increasing values of the vertical magnetic field strength. For 
magnetorotationally stable equilibria, the potential difference increases faster 
than the fourth power of the magnetic field strength, quickly exceeding a value 
corresponding to the central temperature of the disk, and is minimized with 
respect to i at i f« 38°. We show that this property is relatively insensitive to 
the form of the opacity law. Our results suggest that an additional source of 
energy, such as coronal heating, may be required for the launching of an outfiow 
from a magnetized disk. 

Subject headings: accretion, accretion disks — hydrodynamics — MHD 
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1. Introduction 

The highly coUimated outflows associated with many astrophysical objects, from young 
steUar objects to active galactic nuclei, are widely beheved to originate in accretion disks 
threaded by magnetic fields (see, e.g., Livio 1997 for a recent review). While winds from 
rotating stars such as the Sun are influenced by magnetohydrodynamic (MHD) effects, in 
particular with regard to the magnetic braking of the star (e.g. Mestel 1968), these outflows 
are driven predominantly by thermal effects in the corona of the star. In the case of 
accretion disks, however, the possibility arises of a wind driven predominantly by dynamical 
effects. 

The influential model of Blandford & Payne (1982) established the significance of the 
angle of inclination, i, of the poloidal magnetic field lines to the vertical at the surface 
of the disk. They showed that, when i > 30°, an outflow is driven spontaneously by the 
centrifugal force (as viewed in the frame of reference rotating at the local Keplerian angular 
velocity) . This effect occurs because, for material forced to rotate at the angular velocity of 
the foot-point of the magnetic fleld line, the centrifugal-gravitational potential decreases, 
rather than increases, along the field line leaving the surface of the disk. 

This result should not be interpreted as meaning that a wind can fiow freely and 
steadily from the surface of the disk without any thermal assistance. If that were possible, 
it would lead to the dynamical disruption of the entire disk, and it would be more correct 
to say that the disk was never in equilibrium to begin with. Blandford & Payne (1982) 
appreciated that their model of a 'cold' wind does not properly describe the flow in the 
neighborhood of the surface of the disk, where thermal effects must become important. 
Indeed, the flow must pass through a sonic (strictly, slow magnetosonic) point in this 
region and can do so only with thermal assistance. The purpose of our investigation is to 
determine quantitatively how much assistance is required. 
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We are interested in the internal equilibrium of the disk and in the properties of the 
wind in the region immediately above the surface of the disk. It is therefore sufficient to 
consider a radially localized region of the disk and solve for its vertical structure. The 
present work is most closely related to that of Wardle & Konigl (1993), who apphed a 
similar procedure to weakly ionized, protostellar disks in which ambipolar diffusion is 
important. There are several significant differences in our approach, however, which are 
described in detail at the end of this paper. 

In several other studies (e.g. Cao & Spruit 1994; Lubow, Papaloizou, & Pringle 
1994; Kudoh & Shibata 1997), the shape and angular velocity of the magnetic field lines 
is assumed, but we determine them self-consistently by solving the relevant equations. 
Another approach used recently (e.g. Ouyed & Pudritz 1997) is to solve for the wind flow 
in the corona of the disk using axisymmetric numerical simulations. We emphasize that, 
since these simulations describe the flow only after it has become supersonic, they are 
complementary to, but quite distinct from, the present work. 

Recently, Livio (1997) reviewed the observations of all the classes of astrophysical 
objects with which jets are associated. Working on the assumption that the same 
mechanism of acceleration and coUimation operates in all these systems, he argued that 
the most probable such mechanism is indeed the predominantly centrifugally driven MHD 
wind described above. He also suggested that an additional source of energy is required to 
produce jets that are sufficiently powerful to be observed. In this paper, we will provide 
qualitative and quantitative evidence in support of this hypothesis. 

It is now well established that weakly magnetized accretion disks are subject to a 
magnetorotational instability (Velikhov 1959; Chandrasekhar 1960; Balbus & Hawley 1991) 
of which the non-hnear development is a state of MHD turbulence (Brandenburg et al. 
1995; Hawley, Gammie, & Balbus 1995; Stone et al. 1996). For this reason, we include a 
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stability analysis as part of this investigation. However, we do not attempt to address the 
precise relation between MHD turbulence and MHD winds. 

The remainder of this paper is organized as follows. In §2, we present the equations 
governing the equilibrium of the disk, and describe the numerical solutions in terms of 
'standard', dimensionless units. We also discuss the stability of the equilibria. In §3, we 
consider the dynamics of a transonic outflow in the corona of the disk, and identify the 
potential difference that must be overcome by such an outflow. In §4, we translate the 
results into physical units. Finally, in §5, we discuss the implications of our analysis, with 
an emphasis on cataclysmic variable (CV) systems. 



2. Equilibrium of the disk 

2.1. Basic equations 

In this section we describe the local vertical structure of a magnetized accretion disk 
such as might be found in a CV system. The disk is taken to be a steady, axisymmetric 
MHD flow in the gravitational potential of a spherical mass M. The usual cylindrical polar 
coordinates (r, 0, z) are adopted. In the limit of a geometrically thin disk, the governing 
equations may be simplified considerably by means of an asymptotic analysis (Ogilvie 
1997a; hereafter. Paper 1) in which the small parameter e is a characteristic value of H{r)/r, 
where z = H{r) is the location of the upper surface of the disk at radius r. They then 
take the form of ordinary differential equations (ODEs) in the vertical coordinate z at each 
radius r separately, which describe the local vertical equilibrium of the disk. There is also 
an integral relation describing the global magnetic structure, but here we are concerned 
only with the local aspects of the problem. 

In the present work we do not repeat the formal asymptotic analysis, but present 
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the simplified, approximate equations that result from it. The dynamical aspects of the 
equilibrium were considered in Paper I, but the treatment of thermal and radiative physics 
below is new.0 We consider the equations below to afford a minimal description of the 
problem while allowing quantitative predictions to be made. 



2.1.1. Balance of forces 



The angular velocity of the fluid is written as i7 = f2o + where 

( 



\ ) 

is the Keplerian value, which is independent of z, and Qi is the deviation from Keplerian 
rotation, which depends on z, but is much smaller than Qq. The radial component of the 
momentum equation then becomes 

-2prnon^ = ^^. (2) 
^0 oz 

(We write the permeability of free space as /io, but will substitute /xq = 47r later for 
calculations in CGS units.) This equation states that the radial tension force associated 
with the bending magnetic field produces the deviation from Keplerian rotation. Radial 
gradients of fluid and magnetic pressures are negligible, as are the inertial effects of the 
accretion flow. The vertical component of the momentum equation is 

This equation states that both the pressure of the radial magnetic field and the vertical 
gravitational force act to compress the disk. The pressure of any toroidal magnetic field is 
neglected, as is the self-gravitation of the disk. 



■^The scalings correspond to the 'weakly magnetized disks' of Paper I, although equilibria 
resembling the 'strongly magnetized disks' will appear in the limit of a strong magnetic field. 
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2.1.2. Magnetic field 
The solenoidal condition on the magnetic field imphes 

The induction equation reduces to the condition of isorotation, B-VJl = 0, which becomes 

-^B, + B,fl = 0. (5) 

The efi^ects of any turbulence in the disk on the mean magnetic field are neglected. If the 
scale of the turbulent motions were small compared to the scale on which the magnetic field 
varies, it might be possible to describe the mean turbulent EMF in terms of an a-effect, 
representing the regeneration of the mean magnetic field, and a /3-effect, representing 
turbulent diffusion (e.g. Moffatt 1978). However, there is no such separation of scales in the 
present problem and the effects of any turbulence are completely uncertain. Although these 
efi:ects may be important in practice, our assumption of isorotation is the simplest way to 
proceed in the face of this uncertainty. 



2.1.3. Thermal and radiative physics 

Instead of prescribing a polytropic relation between pressure and density, as was done 
in Paper I, we now include thermal and radiative physics. The disk is assumed to be 
optically thick, with Rosseland mean opacity n. The vertical radiative energy flux is then 

16aT^ dT 



3kp dz ' 

and the energy equation may be written 



(6) 



dF 3^ , , 

— _ -iZoCT^^, [7) 
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where is the component of the stress tensor responsible for angular momentum transport 
and dissipation of energy. (Note that the solutions we are considering are stable against 
convection.) 

For a complete specification of the problem, an equation of state must be supplied, and 
also prescriptions for the opacity and stress. We use the ideal gas law, 

kpT 

P = , 8 

where /x is the mean molecular weight of the gas. This assumes that the radiation pressure 
is negligible. For the stress, we adopt the conventional a-prescription, 

ar4> = ap, (9) 

where a is a constant. For the opacity, we assume a Kramers law, 

K = n^pT-"\ (10) 

where kq is a constant. This is an adequate approximation when free- free absorption is the 
dominant source of opacity, as is usually the case in CV disks. We consider more general 
opacity laws in §4 below. While equation (^) represents the most uncertain element of this 
treatment, it does allow a comparison to be made with previous work. 

2.1.4- Boundary conditions 
We solve the equations in < z < H, and apply symmetry conditions 

Br. = (11) 

and 

F = (12) 
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at 2; = 0. The magnetic field then has dipolar symmetry. The boundary conditions at 
z — H are taken to be 

p = 0, (13) 

T = (14) 

and 

Bj. — S^tani, (15) 

where i is the angle of inclination of the magnetic field to the vertical at the surface of the 
disk. 



2.1.5. Conservation of angular momentum 

A further relation, which allows quantities to be expressed in terms of the accretion rate 
M, is obtained from the conservation of angular momentum. Making the usual assumptions 
about the nature of the boundary layer near the surface of the central object at radius r* 
(e.g. Pringle 1981), we obtain 

r^QofM = 2Tir^ J cr^^ dz, (16) 



where 



Then 



/ = i-(7) ■ (17) 



= ^^IfM, (18) 



where the subscript 's' denotes the value at the surface of the disk. 

This is one equation in which a toroidal magnetic field might have a significant effect. 
If the disk has an outflow which exerts a torque comparable to that due to cr^^, it can be 
considered to increase the effective value of Fg in this equation. However, it is assumed that 
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the rate at which mass is lost to the wind is small compared to M. (See, e.g., Livio 1997 
for observational evidence suggesting that this is indeed the case.) 

2.2. Transformation of the equations 

We now recast the equations in dimensionless form by means of the following 
transformations . 

z ^ zH, (19) 
= {^^ Qo, (20) 



T = TH^nlO^'^, (24) 



The dimensionless equations are 



(26) 

-pz + — ~ , (27) 



dz 2B 
dp 

dz B 



z 



dBr _ 2pQi 
dz ~ B, 



(28) 



^ = -p^f-'^/^F, (29) 
dz 

dP 

V = PT. (31) 



We then have a fifth-order system of non-hnear ODEs on < z < 1, with boundary 
conditions 



4(0) = 0, (32) 

F(0) = 0, (33) 

P(l) = 0, (34) 

T(l) = 0, (35) 

Br{l) = S^tani, (36) 



and with two dimensionless parameters, Bz and i, which may be taken to be non-negative 
without loss of generahty. Note that the dimensionless viscosity parameter a has been 
scaled out of the equations. We remark that a similar reduction of the equations could be 
made for any equation of state, stress prescription and opacity law that are 'simple', in the 
sense of being monomials of the thermodynamic variables. 

Dimensionless variables written with tildes will be said to be expressed in 'standard' 
units. These units allow the simplest and most natural presentation of the equations and 
their solutions. They are based on the length scale H, the time scale Q~^, and mass and 
temperature scales derived from the coefficients appearing in the equation of state and the 
opacity law; they also incorporate the scalings of the solutions with H/r and a. However, 
since H is not known a priori, a further transformation is made in §4 below to express 
quantities in 'physical' units based on the accretion rate. For this reason, the interpretation 
of numerical results in the next section should be regarded as provisional. The final 
interpretation is offered in §4. 
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2.3. Numerical solutions 

An expansion of the required solution of equations p6|)-(|3T[) in powers of (1 — 5) is 
readily obtained which allows the equations to be integrated smoothly out of the singular 
point z = 1 towards z = 0. Shooting is then required in two dimensions to match the 
symmetry conditions on z = 0. As discussed by Ogilvie (1997b), for certain values of the 
parameters there exist solutions in which the magnetic field bends more than once as it 
passes through the disk. Such 'irregular' equilibria are known to be unstable and are not 
discussed here. The regular equilibria occupy a connected region in the parameter space, 
as shown in Figure 1. The edge of this region has a form similar to that obtained for 
polytropic disks. Equilibria with any angle i < 30° exist for any value of Bz, but angles 
i > 30° can be achieved only if the magnetic field is sufficiently strong. As the edge of 
the solution manifold is approached, the effective gravitational acceleration parallel to the 
magnetic field at the surface of the disk tends to zero, which in a time-dependent situation 
would lead to the dynamical disruption of the disk. 

EDITOR: PLACE FIGURE |T| HERE. 

The equilibria exhibit a wide variety of behavior in different parts of the parameter 
space, and certain features of this should be explained. First, when i = 0, the equilibria 
are unaffected by the magnetic field and all take the form of the unmagnetized solution. 
Secondly, when i is fixed, with < i < 30°, and the limit i?^ ^ is taken, the equilibria 
eventually become physically unrealistic. Much of the disk becomes almost evacuated, 
with mass concentrated near the equatorial plane and near the surface. This behavior is 
the only way to prevent the magnetic field lines from bending many times as they pass 
through the disk, when the field is very weak. While these equilibria are almost certainly 
unstable to overturning, this behavior occurs only when they are already unstable to the 
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magnetorotational instability, as described in §2.4 below. Thirdly, when i is fixed, with 
i > 0, and the limit ^ oo is taken, the equilibria approach a strongly magnetized limit 
in which they are compressed only by the Lorentz force and not by gravity. We find the 
asymptotic form of the solution to be 

= -Bl''/'''{tani)-y'''X + 0{l), (37) 

p{z) ~ (5,tan026A9yi(5), (38) 

p{z) ~ {B,tanify2{z), (39) 

Br{z) ~ {B,ta,ni)ys{z), (40) 

f{z) ~ iB,tamy'/''y,{z), (41) 

F{z) ~ (B.Unify^iz), (42) 

where A is an eigenvalue of the non-linear system 

y'^ = -2Xym, (43) 

y's = 2A2/1, (44) 

v'a = -vlvl^^'^Vb, (45) 

= y2, (46) 

y2 = yiy4, (47) 

with boundary conditions 

^3(0) = 0, (48) 

VM = 0, (49) 

yi(l) = 0, (50) 

y4(l) = 0, (51) 

1/3(1) = 1. (52) 
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Here, a prime denotes differentiation with respect to the argument. The numerically 
determined solution is A 1.916. Also required is the quantity 1/5(1) ~ 0.1607. This 
limiting behavior of the equilibria is used in §4 below. 

The plasma beta, being the ratio of the gas pressure to the magnetic pressure, is often 
taken to be a dimensionless, inverse measure of the strength of the magnetic field. However, 
as explained in Paper I, this is not appropriate in the case of strongly magnetized accretion 
disks. In the limit S^. — > 00, the plasma beta, evaluated on the equatorial plane, tends to a 
constant value which depends on i. This is because the nature of the vertical equilibrium in 
the strongly magnetized limit requires that the gas pressure increase in proportion to the 
magnetic pressure. 



If the magnetic field is sufficiently weak, the equilibria are expected to be unstable 
to the magnctorotational instability (Balbus & Hawley 1991). The instability of a thin 
disk containing a bending poloidal magnetic field has been analyzed by Ogilvie (1997b), 
who found that the curve of marginal stability in the parameter space could be located by 
solving the equations for an equilibrium possessing a mode with zero frequency and zero 
radial wavenumber. In terms of the Lagrangian displacement ^, these are 



2.4. Magnetorotational stability of the equilibria 




(53) 



and 




dSU 



o2 ^'^^ 



(54) 



dz 



where 




(55) 
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is the Eulerian perturbation of total pressure. This assumes that the perturbations are 
adiabatic, with adiabatic exponent 7. The relevant mode is the first mode of odd symmetry, 
which satisfies the boundary conditions 

= ^ = (56) 



at z = 0, and 



^-tan.^ = (57) 



at z = H. The dimensionless form of these equations is identical except for the inclusion of 
the tildes and the omission of f^o- 

The curve of marginal stability, for the case 7 = 5/3, is shown in Figure 2. Again, this 
is qualitatively similar to curves obtained for polytropic disks. 



EDITOR: PLACE FIGURE | HERE. 



As discussed by Ogilvie (1997b), it is expected that the equilibria that are stable 
to the magnetorotational instability are also stable to the magnetoconvective (Parker), 
interchange, and bending instabilities, unless the magnetic field provides most of the 
support against the radial gravitational force. Although a weak, global, non-axisymmetric 
instability may remain, this is unlikely to be dynamically important in a thin disk. 



3. Dynamics of an outflow in the corona of the disk 

Following Paper I, we now analyze the dynamics of a transonic outflow in the region 
immediately above the disk, which we refer to as the 'corona', and which is defined by 
H < z <^ r. The density of the wind is very much smaller than that of the disk, possibly by 
a factor 0(e"^), and to match the solutions in detail would require high-order asymptotics 
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well beyond the scope of this analysis. Instead, we consider that the disk acts as a reservoir 
which can supply any reasonable mass flux to the wind, the corresponding velocity field 
required in the disk being extremely small. 

In the corona, the magnetic field is force-free to a very good approximation, and the 
field lines are straight on the length scale H. They act as rigid channels for the wind and 
also enforce isorotation. The dynamics of the outflow in this region depends critically on 
the centrifugal-gravitational potential which may be computed as follows. Since the 
magnetic field is force-free, dB^jdz = 0, so that Bj. = Bztani throughout the corona. 
Equation (^) may then be integrated to give 

ill = ihs H , (58) 

2r 

where, again, the subscript 's' denotes the value at the surface of the disk. The effective 
gravitational acceleration is 

g = 2ri7of^i — ^Iz e^, (59) 

and its component measured parallel to the magnetic field (and towards the surface of the 
disk) is 

g\\ = —(Standi — l)Ql{z — H) cosi + Qq^QqH — 2i7isrtanz) cosi. (60) 

If, as we assume, i > 30°, then g\\ decreases linearly with increasing z and goes to zero at 
z = z^onic, given by 



2;sonic = H ^ • TT7^ • (61) 



{floH — 2r2isrtani 
(Standi - 1)^0 

This name is appropriate because, as shown below, z = -Zsonic is the expected location 
of the sonic point of a transonic wind. Consider a single magnetic field line, and let z 
parametrize the position along it. When gfy is integrated along the field line, it yields the 
centrifugal-gravitational potential 

= _i(3tan2z - l)Ql{z - ^sonic)', (62) 
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as measured from the sonic point. In particular, the potential difference between the sonic 
point and the surface of the disk is 



{noH - 2^]lsrtan^)2 
2(3tan2i - 1) 



A$ = ' Z~ 2 ;r • (63) 



When i approaches 30° (from above), the height of the sonic point increases without 
bound, as does the potential difference. However, these expressions are valid only if the 
height of the sonic point calculated from equation (^) satisfies 2;sonic r, since otherwise 
the quadratic approximation to the potential ceases to be valid and the curvature of the 
magnetic field lines should also be taken into account. Within this approximation, the 
potential difference is always small compared to GM/r. (By comparison, when i < 30°, the 
potential difference is comparable to GM/r.) 

We now consider the dynamics of a transonic wind flowing along the magnetic field 
lines. The wind is treated as isothermal, since the optical depth is presumably small in the 
corona, and therefore 

p = c'p, (64) 

where 

(65) 



is the isothermal sound speed. Mass conservation requires that the mass flux density 

pu = J-" = constant (66) 
be constant along the flow. The Bernoulli equation states that 

+ In p + ^'^s = constant (67) 
is also constant following the flow. On differentiating this equation we obtain 

(^2 _ c^)^^ = (3tan2z - l)nl{z - z.onic), (68) 
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which demonstrates that the sonic transition must occur at z = Zgomc- In terms of the Mach 
number Ai = u/c, we then have 

1{M^ -l)-lnM = -$=Vc2, (69) 

and, in particular, 

liMl -l)-lnM, = A$/c2 (70) 

is a transcendental equation for the Mach number A4s at the surface of the disk. The mass 
flux density is then 

^ = MsPsC, (71) 

where ps is the density of the wind at the surface of the disk. When A$ ^ c^, an 
approximate solution is 

>ls~exp(-A<l>/c2-i), (72) 

which implies 

J^^p,cexp{-A^/c^-l). (73) 

The numerical solution of equation ( |70D in Figure 3 shows that this approximation is 
accurate even for moderate values of A$/c^. 

EDITOR: PLACE FIGURE | HERE. 

Evidently the meaning of equation (^) is that the outflow is severely choked if the 
potential difference is much larger than c^. For any equilibrium we can compute A$ and 
deduce the sound speed - and therefore the temperature - required in the corona for the 
outflow not to be suppressed. The height of the sonic point and the potential difference are 
plotted as functions of i for various values of in Figure 4. These equilibria are all stable 
to the magnetorotational instability. As i increases from 30°, the sonic point approaches 
the surface of the disk and the potential difference decreases, until the equihbrium ceases to 
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exist. When the magnetic field is stronger, the sonic point is more distant and the potential 
difference is larger. However, for any field strength, an angle i can be found for which the 
potential difference is arbitrarily small. In the next section we show that this property 
is lost when the solutions are expressed in more physical units which take account of the 
compression of the disk by the magnetic field. 

EDITOR: PLACE FIGURE | HERE. 



4. The solutions expressed in physical units 

4.1. Definition of units 

Although the dimensionless forms defined in §2.2 are convenient for the purposes of 
computation, one would like to express the magnetic field, for example, in units of gauss 
rather than in the dimensionless form B. The relation between B and B involves H which 
is not known a priori. However, the accretion rate M is determined a priori in the sense of 
being of a global property of the disk, independent of radius, whose value can be estimated 
observationally. From the numerical solutions, we know Fg for each equilibrium, which 
allows us to write 



combining equations pq ) and (^). We then have 



where 

U, 
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is a suitable unit of length. Similarly, we may write 

B = BF-i^/^° Ub (77) 



and 



where 



T = rF-3/^° Ut, (78) 



and 



..^(^)'""(f)''"^«r(^)-(S^)"'"°/— ,so) 

are suitable units of magnetic field strength and of temperature. 

In terms of CGS units, we find 

Uh ^ 1.2xlO«a-ViOMr'/'i??^V'/'°Mfrcm, (81) 
Ub ~ 1.0 xl03a-9/20Mi'/^'i?ro''^'V''/'°MiY^'° gauss, (82) 
Ut ~ 1.3xl0^a-V5My'i?ro'V/'°Mf/'°K, (83) 

where Mi = M/Mq, Rw = it:/(10^° cm), Mie = M/(10^^ gs"^), and we have used the 
values iJ, ^ 0.6 and Kq 6.4 x 10^^ cm^ g~^ (cf. Novikov & Thorne 1973) appropriate for CV 
disks. 



4.2. Solutions with a purely vertical magnetic field 

When the magnetic field is purely vertical, it does not affect the equilibrium. The 
solution is the same for all values of Bz, and is found numerically to have Fg ^ 0.0007041. 
We then obtain 

H ^ 2.971 Uh ^ 3.5 x 10« a'^/^'^ M^'^' R%' f/^^ mU'' cm. (84) 
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We note that this is larger by a factor of approximately 2 than the value quoted by, for 
example, Frank, King, & Raine (1985), which was based on order-of-magnitude estimates 
(although we define H to be the true semi-thickness rather than an approximate scale 
height). The magnetic field strength for marginal magnetorotational stability is found to be 
^ 0.06327, or 

^ 1.383 Ub ^ 1.4 X 10^ «-9/20^7/16^-^21/16^17/40^17/40 ^^^^^^ ^g^^ 

If the magnetic field results from the aligned dipole field of a central white dwarf of radius 
5.0 X 10^ cm, a value of 1.4 x 10^ gauss at a radius of 10^° cm in the disk would correspond 
to a value of 2.2 x 10'^ gauss at the poles of the white dwarf. This means that the disk 
would be unstable at this radius in many cases of astrophysical interest, but may be stable 
in some systems such as V1500 Cyg. 

4.3. Solutions with an inclined magnetic field 

We now redraw Figure 2 in physical units, using BzF~^'^^'^^ as the ordinate rather than 
Bz- Unfortunately the resulting graph (Figure 5), while more physically meaningful, is more 
difficult to interpret. The mapping from B^ to B-^F"^'^/'^^ is not one-to-one, with the result 
that the solution manifold folds over on itself on the left-hand side of Figure 5. However, the 
'folded' solutions are unstable and this detail need not be pursued here. More interesting is 
the way in which the curve on which the equilibria cease to exist is transformed relative to 
Figure 2. This distortion occurs because equilibria with highly inclined magnetic field lines 
are strongly compressed by the Lorentz force and therefore have a higher pressure than 
the unmagnetized solution with the same H . The stress and torque are correspondingly 
increased, so that is large, and therefore B^/Ub is smaller than might be expected. 



EDITOR: PLACE FIGURE | HERE. 
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The existence of the fold means that, when the magnetic field strength is reduced 
below the stability boundary, the equilibria eventually cease to exist. However, it is of little 
importance whether, in this part of the parameter space, no stable solution exists or no 
solution whatever exists. The result is likely to be a state of MHD turbulence. 

The height of the sonic point and the potential difference, expressed in physical units, 
are plotted as functions of i for various values of Bz/Ub in Figure 6. The height of the sonic 
point is compared with if, both expressed in units of Uh- For the potential difference we 
introduce a physical unit 

[/^ = ^ 1.8 X 10^2 a-^''Ml"R^^"fl''MU''' cm' s'^. (86) 

Then A$/f/$ can be compared with Tc/Ut, where Tc is the central temperature of the disk. 

EDITOR: PLACE FIGURE | HERE. 

In the range 1.55 <^ Bz/Ub <^ 1.93, there exist stable equilibria such that, when i is 
increased from 30°, the potential difference falls rapidly from infinity to zero as the sonic 
point approaches the surface of the disk, and then the equilibria cease to exist. This 
behavior is similar to the interpretation offered earlier, based on Figure 4 and standard 
units. There are also equilibria with more highly inclined magnetic field lines. However, 
for Bz/Ub ^ 1.93, a different behavior is found. As i is increased from 30°, the equilibria 
continue to exist for all angles up to 90°. The potential difference has a minimum, typically 
in the range 38° <^ i <^ 43°, and then increases again. The sonic point at first approaches the 
surface of the disk, but then recedes from it. Moreover, as Bz/Ub is increased, the minimum 
potential difference increases very rapidly and quickly exceeds a value corresponding to the 
central temperature of the disk. 
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4.4. Limiting behavior for strongly magnetized disks 

When Bz/Ub is sufficiently large, the limiting form of the equilibria given in §2.3 
applies. We then find that the potential difference is 

^\Ub) (3tan2z-l)' ^ > 

where 

Ci = 2A2[|/5(1)]30/i9 ^ 0.4095. (88) 

This increases very rapidly with increasing magnetic field strength. The minimum with 

respect to i occurs at 

/l4 

i = arctan W ^ « 37.96°. (89) 
Similarly the height of the sonic point is 

cfj.y^^" /^--r". (90) 



Uh \UbJ (3tan2i-l)' 

where 

C2 = 2A[2/5(1)]^^/^^ ~ 0.9050, (91) 
and this has a minimum with respect to i at 

i = arctan ^ 74.05°. (92) 

To show that these power laws and characteristic angles are relatively insensitive to the 
opacity law, the analysis can be repeated using a more general relation of the form 

K = «op^T^. (93) 

We then find that the potential difference is 

A$ /B^\ 4(5+23;-2/)/(5+^-2/) (tan 2)4(5+2x-:y)/(5+x-2/) 

'^i 771 T^TTZX^ — ' (94) 



\UbJ (3tan2i-l) 
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where Ci depends on x and y, and of course the units f/$, Ub, etc., are differently 
but analogously defined. This shows that the potential difference typically scales as 
approximately the fourth power of the magnetic field strength. The minimum with respect 

to i occurs at 

i = arctan , 



2(5 + 2.-,) 



provided that 

5 + 3x — y 



5 + X — y 

Similarly, the height of the sonic point is 



\\ 3(5 + 3x -y)' 

> 0. (96) 



J^_^s2i5+2x-y)/i5+x~y) ^_^^^^y{5+2x-y)/{5+x-y) 

•^2 ( 77^ /o.„„2 ■ T{ ' (97) 



Uh \UbJ (3tan2i-l) 

and this has a minimum with respect to i at 



provided that 



5 + 2x-y 

arctan A / , [9t 

V 3x 



^ > 0. (99) 



5 + X — y 

Criterion ( p6D is satisfied for most reasonable opacity laws, as shown in Figure 7. In the 
case of electron-scattering opacity, for which x = y = 0, the potential difference scales as 
{Bz/UbY is minimized with respect to i at i ~ 39.23°. Criterion (p9|) marginally fails 
to be satisfied, however, so that the height of the sonic point continues to decrease as i 
approaches 90°. 

EDITOR: PLACE FIGURE HERE. 

In Figure 8 we compare the true minimum of the potential difference with respect 
to i with the value predicted from equation ([87|). This shows that the limiting behavior 
described in this section is achieved very rapidly as the vertical magnetic field strength is 
increased beyond the stability boundary. 
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EDITOR: PLACE FIGURE | HERE. 
5. Discussion 

In this paper we have solved for the local vertical structure of a magnetized accretion 
disk such as might be found in a CV system. The magnetic field was assumed to enforce 
isorotation, and the deviation from Keplerian rotation was taken fully into account. Angular 
momentum transport in the disk was parametrized by an a-prescription in the conventional 
way, and a Kramers opacity law was assumed to hold. We have shown that, when quantities 
are expressed in physical units, equilibria that are magnetorotationally stable can be found 
in which the angle of inclination, i, of the magnetic field to the vertical at the surface of the 
disk has any value in the range < i < 90° if the magnetic field is sufficiently strong. 

We have analyzed the dynamics of a transonic outflow in the corona of the disk when 
i > 30°, and, in particular, have shown that a certain potential difference must be overcome 
by such an outflow. When the equilibria are very close to the magnetorotational stability 
boundary, the potential difference is relatively small and can in fact be made arbitrarily 
small by approaching the edge of the solution manifold. For more strongly magnetized disks, 
however, the potential difference increases faster than the fourth power of the magnetic field 
strength, and is minimized with respect to i at i ^ 38°. These properties are relatively 
insensitive to the opacity law. 

We have used a local representation which neglects the global form of the disk and 
wind, focusing instead on the vertical structure at a single radius. In this respect our 
analysis is comparable with that of Wardle & Konigl (1993). However, in the limit in which 
ambipolar diffusion is negligible, as is the case in the accretion disks of many classes of 
objects such as CVs, no meaningful solution of their equations can be obtained, even though 
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some of their equations correspond in this hmit to the ones we have solved. To insist that 
the disk achieves a balance between the outward diffusion of magnetic flux and its inward 
advection by the accretion flow is certainly attractive, but constrains the problem in such a 
way that it appears that diffusion is driving the outflow. 

Certain other features of our analysis should be emphasized. First, we do not assume 
that the wind accounts for all, or even most, of the angular momentum transport in the 
disk. In the case of CV disks, for example, we can be confident that the agent supplying 
most of the torque on the disk is not an MHD wind; see the discussion by Livio (1997). 
Secondly, we consider the disk to have a well-defined surface which is hardly affected by 
the presence or absence of an outflow above it. We do not attempt to make a smooth 
transition between the disk and the wind, except in requiring the continuity of the mass 
and magnetic fluxes, because the density of the wind is so small that a detailed matching 
procedure would require high-order asymptotics. Thirdly, we do not consider ambipolar 
diffusion, Ohmic resistivity or turbulent diffusion of the magnetic field. Instead, we assume 
that isorotation holds on magnetic field lines and we allow for a slow accretion of magnetic 
flux by the disk. Finally, we have considered an optically thick disk rather than assuming 
the disk to be isothermal. We note that, if the disk is isothermal, it is easier to obtain an 
outflow of reasonable strength without postulating a hot corona, because the region above 
the disk always has the same temperature as the central temperature. 

The interpretation of our results depends to some extent on the class of accretion disks 
under consideration, and on the source of the mean magnetic field. In the case of CV disks, 
for example, the magnetic field of the white dwarf can be important and may be sufficient 



to make at least the inner part of the disk magnetorotationally stable (cf. eq. [Q). In that 
case, even if the magnetic field lines are inclined such that i > 30°, our results suggest 
that the potential difference may be too large to allow a significant outflow unless the 
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disk has a hot corona or access to an additional source of energy, such as nuclear burning, 
in accordance with the hypothesis of Livio (1997). However, the interaction between a 
magnetized central object and the disk may be more complicated than we have allowed for 
in this paper (cf. Livio & Pringle 1992; Miller & Stone 1997). In other cases, magnetic flux 
may be accreted from the environment or may be a remnant of the formation of the disk. 
In disks that arc magnctorotationally unstable, the dynamics of the mean magnetic field is 
much less certain, but it may be possible for these systems to regulate the distribution of 
magnetic flux so as to remain close to the stability limit, and thereby avoid incurring a very 
large potential difference. 

GIO thanks the Space Telescope Science Institute for its hospitality. ML acknowledges 
support from NASA Grant NAGW-2678. 
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Fig. 1. — Parameter space of equilibrium solutions. Bz is the vertical magnetic field in the 
disk, expressed in standard units (eq. p3|), and i is the angle of inclination of the magnetic 
field to the vertical at the surface of the disk. Dashed line: edge of the solution manifold; 
regular equilibria exist everywhere to the left of this line. Dotted line: drawn at an angle 
i = 30°, which is critical to the analysis in §3. (In this graph only, a logarithmic scale of 
magnetic field strength is used in order to give a broad perspective on the parameter space.) 

Fig. 2. — Solid line: stability boundary to the magnetorotational instability, for the case 
7 = 5/3. Dashed line: edge of the solution manifold, as in Figure 1. (Note that the axes are 
quite different from those of Figure 1.) 

Fig. 3. — Dependence of the Mach number Ais of the wind at the surface of the disk on the 
potential difference A$ between the sonic point and the surface, expressed in units of c^, 
where c is the isothermal sound speed. The solid line denotes the true solution, while the 
dashed line represents the approximation given in equation (|72D. 

Fig. 4. — Variation with i of the height of the sonic point (a) and the potential difference 
between the sonic point and the surface (b), both expressed in standard units, for equilibria 
with Bz = 0.1 {solid line), 0.2 {dotted line) and 0.3 {dashed line). 

Fig. 5. — Parameter space with the vertical magnetic field Bz expressed in physical units (eq. 
|79|). Solid line: stability boundary. Dashed line: edge of the solution manifold. Dot-dashed 
line: fold in the solution manifold. No regular equilibria are found below a curve which 
consists of the dot-dashed hne continued by the dashed line. 

Fig. 6. — Panels a, c, and e: height of the sonic point {solid line) compared with H {dotted 
line), both expressed in physical units, as functions of i, for equilibria with Bz/Ub = 1.8 
(a), 2.0 (c), and 3.0 (e). Panels b, d, and f: potential difference {solid line) compared with 
the central temperature {dotted line), both expressed in physical units, for equilibria with 



Bz/Ub = 1-8 (6), 2.0 {d), and 3.0 (/). Note the different scale for the potential difference 
and temperature in panel /. 

Fig. 7. — Parameter space of opacity laws k = Kop^T^. The points marked 'K' and 'ES' 
correspond to Kramers opacity and electron-scattering opacity, respectively. Solid lines: 
contours of the angle of inclination for which the potential difference is minimized when 
expressed in physical units (eq. [0). The contour values are 31°, 32°, . . . , 40°. Dotted lines: 
contours equivalent to 30° and 90°; in the narrow sector between these lines, a minimum is 
not obtained. 

Fig. 8. — Solid line: true minimum of the potential difference with respect to i, as a function 
of the vertical magnetic field strength, both expressed in physical units. Dotted line: the 
limiting form calculated from equation (^). 
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